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Suiamary: The investigations - in an earlier publication - con- 

cerning force-deflection diagrams of airplane shock-absorber struts with 
linear damping are repeated with the assumption that the damping be pro- 
portional to the square of the spring-compression velocity. The dia- 
grams with quadratic damping essentially differ from those with linear 
damping. The comparison with drop-hammer diagrams that had been plotted 
shows that the damping actually depends linearly, or at least almost 
linearly, on the spring-con^iression velocity. 

The unsolved problems in \mdercarriage construction are still so 
n\mierous that, like a hydra, from each problem that has been solved, 
two new ones grow out. Thus my present investigation is an outgrowth 
of a previous report (ref. 1) on the comparison of "bviffered" shocks 
with "mbviffered shocks"; however, it encompasses only a single addi- 
tional problem. On the other questions that arose I hope to report in 
sequence later on.l 

The present investigation thus repeats the considerations reported 
in said treatise with the difference that here a damping is assumed 
which is proportional to the square of the spring compression velocity, 
while there we dealt with linear damping. While I have found no direct 
reference in the literature as regards to the type of damping expectable 
in airplane spring struts, personal experiences and oral communications 
from outside sources would indicate that the linear damping was set up 
mainly with respect to mathematical treatment and only as a second 
thought with respect to physical correctness. Thus the calculatory 
treatment of double-springings, such as spring struts with quadratic 
damping and tires, offers insurmountable difficulties, whereas with 
linear damping an exact mathematical treatment is possible (ref. 2) . 


*"Zur Kenntnis der Kraftwegdlagramme von Flugzeugfederbeinen, 

1. Teilbericht; Vergleich von Diagrammen mit linearer und quadratischer 
Dampfung." Technlsche Berichte, Bd. 11, Heft 2, Feb. 15, 19^4, pp. 51-55. 

^NACA editor ' s note ; Subsequent papers in this series are included 
in this translation as Second and Third Partial Reports. 
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In the following, therefore, the force-deflection diagram with q\iad- 
ratic damping shall be compared with linear damping in order to gain a 
perspective of the physical relationships therein. 

The energy equation for the first spring shock (using the same 
assumptions and designations of the previous report) is 


J V df =J(p^ + cf + k 2 V^^df 

= ^(^0^ - - L)f 

From eq\iation (l) with 




c = ^ = 2o£^ 

M' Mg 


Mg - L - Py _ 


Mg 


— 1 — — cx»^ ^ 


V = (p % 

cu 




one obtains, after differentiation with respect to f 

^ + !!^ v2 + 2«^f - 


df 


2ga = 0 


( 1 ) 


( 2 ) 


(5a) 


✓ 




or 


+ 2(^|» - a) = 0 (5b) 

d\|( 

The general solution (ref. 5) of the differential equation (51^) 

reads 


cp2 = Ae-^i + Bt + C 


(M 
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When inserting the solution (4) into the initial equation (5h) we 
obtain 


Ae-5'l' (- 5 + 4 ^ 2 ) + \|»(4^2® + + (b + - 2aj = 0 


( 5 ) 


The parentheses must disappear if equation (5) should be fulfilled for 
every value of >lr. Thus is obtained 


6 = 4i92 


B = . 


2*62 


C = (a + ^ - 

21^2 \ ^*'^2 




The constant A is derived from the initial condition 




as 


( 6 ) 


(7) 


A = 


/« 4 . 1 

• Ot. + T— — 


form 


^0 - • 4,2^ 

With equations (6) and (8), the solution (4) is written in final 


( 8 ) 




_i- i|( + -i-[a + T-i 
2^0 2^0 V 4d 


( 9 ) 


When setting approximately 
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which is permissible for low dampings, we obtain 


» <Po= - 5^ ♦ 


( 11 ) 


and with 


finally 



( 12 ) 


(15) 


Thus, for not too great dampings, the square of the spring compres- 
sion velocity decreases linearly with the spring stroke . 

2 

In order to determine the maximum deflection \lfg, cp = 0 is 
inserted into equation (9)^ and thus obtain 



( 11 ) 


an equation which may be solved by trial only. 


The maximum spring stroke may be determined with another method of 
approach as well. 


J = s dt = ^ + )l(g(l - 

J Q ^ 


A) 


When calculating the integral, we obtain 


( 15 ) 


J = 





2^2<Po^ - 







(16) 
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According to equation (2) 


0 + a = 1 - A 


and therewith one obtains, eqiiating equations (15) and (16), 
equation (l4) . 

As already established in equation (16), the force S 


P Pv + cf + kpv2 

'Mi* = c + 2.2^2 


The maximum value Sg is obtained by a zero-setting of the 
quotient. We then have 


d\|f 




= 1 + = 0 


or 


\ d\|r /* 2’d2 


When inserting equation (20) into eqiiation (3b), we obtain 




" ' 4(4 - (♦* - “)) 

The insertion of equation (21) into equation (4) leads to 


4(4 ■ ^ 

or, according to equation (6), in general 


B\|r 


+ C = -i-U 

2^2 \H 


- (i|r - a)j 


(IT) 

likewise 

is 

(18) 

differential 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 


(25) 
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to the condition 


= 0 (24) 

This implies that S steadily increases further with \|r and finally 
reaches its optimum for \|f = ». For that reason the highest value of 
S is expectable in the force -deflection diagram at \|f = or f = fg, 

respectively. Thus we have 


Sg — cf + i|(g (25) 

Finally, the springing effectiveness tj is given by the equation 


n = 


2t 


+ 1 - A 


g 


a + 


(26) 


The component of the force S caused by damping is, according to 
equation (I8) 


D = 2^2^?^ = Pq 




2Do, 




2Do, 


(27) 


when the damping force at the start of the spring shock is expressed 
thus : 


Dq = 2i82Cpo^ (28) 

From equation (27) > a = 0, that is, A + a = 1, a variation of the 
damping force D is obtained as shown in figure 1. It is seen that, 
for finite values of the damping also, the course of D as function 
of \|r does not depart essentially from a straight line . With the 
assumption that 


Do - 


a + 



= 0 


(29) 
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or 


Do - Docrit - ^ 



(50) 


D as a function of the spring stroke is exactly a straight line for 
which the equation reads 


D = Docrit - + 


(51) 


For a = 0 this occurs for DQcrit “ shown by the line for 

Do = 5 and cpQ^ = 50. When Dq < Docrit# damping cvirve is concave 
toward the coordinate origin ( Dq = 1 and 2 , cpo^ = 10 and 50). If 
Do ^ Docrit^ "then the ciirve is bent toward the zero-point (Dq = 5^ 
cpo2 = 10> 


From equations (l4) and (28) we obtain, with cpQ^ = 50 as well as 
with 0 = 0 and X = 1^ that is, a = 0, the variation of i^g as a 
function of the initial force Dq, as shown in figure 2. The, latter 
serves as a foiindation for plotting the force-deflection diagrams of 
figiire 3* 


Corresponding to the linear damping expectable with airplane spring 
struts, we obtain with <Pq 2 = 50 and = O.25, from my previous report 

(ref. 1, eqs. (l3a), (l7a), and (20)) 


\|/gl = 5.0515 
% = ^t.525 

Sgi = 5.7525 


(52) 


The diagram d plotted for linear damping was made with the above 
values . As a comparison therewith were shown the three disigrams a, b, 
and c for quadratic damping. 

Diagram a shows the same maxiimm force Sg as the diagram d, 
whereas diagram b begins with the same damping force Dq. Finally, dia- 
gram c has the same maximum spring stroke as diagram d. In figvire 2, 
the points defining the relationship between ijrg and Dq are indexed 

with the letters of the diagram. 



1-8 


NACA TM 1575 


T’.'US; the form of the force-deflection diagrams with quadratic 
damping is essentially different from those of linear damping. The dia- 
grams obtained in drop tests, however, have a similar appearance as dia- 
gram d, from which we may conclude that the actual damping of airplane 
spring struts is exactly, or at least very nearly, proportional to the 
spring compression velocity. 

According to this understanding the quadratic damping thus has 
probably not too great a practical importance and it has been found suffi- 
cient to calculate further on with the compression spring velocity course 
as was given in equation (15) . 


For unchanged energy absorption then the approximation value 
is derived from the equation 


2d + \^gn + Dq , 

P ^gn 



as follows 


(53) 


ign - - 




(5i^) 


The error, introduced by calculating with equation (34) instead of the 
exact equation (l4), is illustrated in figure 4. Here we see that for 
mild damping (E^ < 5 . 5 ) and for cpo > 3 corresponding to the system 
constants of larger aircraft (ref. 1, p. 151 ); the error remains below 
4 percent. For this reason we may calculate from here on without any 
objections with the approximation formula (54) . 


The characteristic chart of the following diagrams with linear 
damping are taken from the oft-mentioned prior report (ref. l) , whereby 
= 0.2 was chosen. Identical energy absorption assumes, according 
to equation (15)> with equal shock velocity Vq also equal maximum 

spring strokes, so long as A 7 ^ 1 , that is, as long as the weight cancel- 
lation (by lift) is incomplete . With full weight cancellation (A = l) , 
however, various comparison diagrams may be foimd with quadratic damping, 
as shown in figure 3 - "the following, with full weight cancellation, 
identical maximum forces are presupposed. 
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When establishing that 


' Wgi’ ^g ' V*gi’ 

^0 = \ 2 p 0 ±> = ^ 2/^1 

wherein the subscript 1 relates to linear damping and subscript 2 to 
quadratic damping, figure 5 is given with a complete lack of weight 
cancellation, whereas figvire 6 is given for fvill weight cancellation. 

The solid line curves relate to nonpreloaded spring struts and the broken 
line curves relate to a pre loading equal to the static load. 

Thus, for^ A = 0, the spring strokes are identical for both damping 
laws, whereas in the practical region of cpQ the maximum forces with 

quadratic damping are 5.5 to 7-5 percent lower than with linear damping. 
The springing effectiveness values ijg are by 6 to 7*5 percent larger 

than the springing effectiveness values t|j^. 


With full weight cancellation (A = 1, fig. 6) the maximum spring 
strokes ilfg2 are, for equal maximum forces, larger than the spring 

strokes by approximately 10 percent; whereas the springing effec- 

tiveness values are, with quadratic damping, on the average 10 percent 
below those with linear damping. 

The initial damping Dq 2 becomes at d = 1 for cpQ = l.i<-8, equal 

to zero (point a) so that the coxorse of the curves for values (p < 1.48 
has no more practical meaning. Thus the QTj-curve ends in point b. 

For smaller values of (Pq, a comparison is no more possible for known 

reasons (ref. l) ; this is, however, of no consequence. 

When, in conclusion, figs. 5 and 6 are viewed again, it may be seen 
that the realization of a quadratic damping law would offer no advantages 
over the linear damping. One should be therefore satisfied with the 
fact that in airplane struts with fluid damping the damping generally 
depends linearly on the spring compression velocity. 

The next partial report, possibly even a third one, shall be devoted 
to the investigation of force-deflection diagrams with a nonlinear 
springing characteristic and linear or quadratic damping. 


Translated by John Perl 
Lockheed Aircraft Corp. 
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Figure 1. 


Figure 2 



Variation of the damping force as function of the spring stroke. 



Do'/fc 

25 *0 


Relationship between maximum spring stroke and initial 
damping (example). 
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Static spring 
characteristic 


Figure 3. - Comparison of three diagrams of quadratic damping (a, b, and c) , 

with a diagram of linear damping (d). 



Figure 4.- Accuracy of the approximation formula 34 compared to the 

exact formula 14. 
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Figures 5 and 6.- Comparison of the maximum forces, greatest spring 
strokes, springing effectivenesses, and initial dampings with quadratic 
damping, with the same values in linear dampii^. 
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SECOND PARTIAL REPORT - APPROXIMATION METHOD FOR THE CALCULATION 
OF FORCE -DEFLECT ION DIAGRAMS WITH A NONLINEAR SPRING 
CHARACTERISTIC AND HNEAR OR QUADRATIC DAMPING* 

By K. Schlaefke 


Sunmary: While undamped nonlinear oscillations on one hand, and 

damped harmonic oscillations on the other may be treated with mathe- 
matical exactness, one is dependent on an approximation method for the 
investigation of damped nonlinear oscillations. In the present report 
will be developed two approximation methods, their use described and 
critically examined. 

Airplane shock-absorber struts, in their most employed form of air- 
spring st.ruts with annular oil damping, represent oscillation elements 
with nonlinear spring characteristic and velocity-proportional damping. 

To the undamped nonlinesir oscillations already a number of thorou^ 
investigations (refs. 1, 2, and 3) have been devoted. With the same 
thoroughness have been explained the damped harmonic oscillations 
(refs. 4, 5, and 6). In contrast thereto nothing as yet has been pub- 
lished on damped nonlinear oscillations that arise out of the combina- 
tion of both above mentioned individual problems. The present investi- 
gation therefore deals with this naturally true oscillation problem. 
Whereby the promise, made at the conclusion of the first partial report 
(ref. 7), shall be fulfilled. 

The calculation of damping-free nonlinear oscillations leads to 
elliptical integrals, whereas that of the damped harmonic oscillations 
leads to differential equations which cannot be solved quite simply. 

With this in both cases the limit is reached for exact mathematical 
solutions. For this reason it is easily understood that a combination 
of both problems remains from the start, not amenable to exact treatment. 
Yet the technical developments require data for construction and test. 

Thus we are* confronted with the problem to find a useful approximation 
method for the calculation of force-deflection diagrams for the first 
landing shock of oil-air spring struts. 

After a few prelimineiry trials I have decided on the method described 
below in which the force-deflection diagreim is subdivided into a greater 


*"Zur Kenntnis der Kraftwegdiagramme von Flugzeugfederbeinen, 

2. Teilbericht: Naherungsverfahren zum Berechnen der Kraftwegdiagramme 

mit nichtlinearer Federkennlinie und linearer Oder quadratischer Dampfung. 
Technische Berichte, Bd. 11, Heft. 4, Apr. 25, 1944, pp. 105-109. 
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number of such small increments that the spring force as well as damping 
force in them may be considered as a straight line with the spring 
compression. Then from the energy-balance an equation may be easily 
derived for the step-by-step calculation of the diagram. 

Figure 1 applies to linear and figure 2 to quadratic damping. The 
designations have been already discussed previously (ref. 5) so that at 
this time we shall, insofar as danger exists in confusing the linear 
with quadratic damping in their appearance side-by-side, designate the 
values with linear damping by the subscript 1, and those of quadratic 
damping with subscript 2. Then the energy balance for linear damping in 
decrement Afj^ is 

i(^-l - "^i^) + Mg(l - X)Afi = Afi 


^i-1 + ^i + ^l^i 


Afi 


( 1 ) 


From this equation is obtained 


Vi = - 


^1 ^i . i/Z^l ^i 


2M 


2M 


) ^ vti - 2g(l - 


( 2 ) 


This is already the formula for the step-by-step calculation of the force 
deflection diagram with an arbitrary spring chart approximated by inter- 
polation and with linear damping. When a step leads to a complex numerical 
value of Vj^, in which case the expression under the root becomes negative, 

calculation is to be made with the formula for the end section. The 
latter reads 


Afn ~ 


M v2 , 
2 n-1 


^n-1 ^n-1 


Mg(l - \) 


A - Jn_i 


Pn-1 + Fn-i 


- Mg(l - X) 


( 3 ) 


2 
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Herein it is approximative ly assumed that the mean force P in the end 
section is equal to the arithmetic mean of Pn_i and this mean 

value is plotted in figures 1 and 2 with a dot and dash line. A desig- 
nates the total work ^ Vq^ and Jn-i denotes the energy absorbed by 
the spring strut up to the spring compression fn_i» 

The energy balance for quadratic damping (fig. 2) is as follows: 


K’^i-l ■ ^i^) + - ^)^i = ^ ^ Af^ 


From this is obtained 


'1-1 -seii- x))Afi 


Vi = 


1 + 


^2 Afj 


The section-formula (3) is also valid for quadratic damping. 


(4) 


(5) 


The application of the formulas shall be shown on an exan®)le. For 
the construction model Bfl09B we have according to Kochanowsky (ref. 8) 

c = 29000 kg/m 1 (6) 

M = 100 kg s^/m J 


With this is derived 


CO = \/| = 17.0294 l/s (7) 

Furthermore we assume (fig. 3 ) 

Fq = 491 kg 
Vq = 4.2 m/s 

= 0.25 

k]^ = 

k^VQ = 3576 kg 

A straight line characteristic chart was purposely assumed so that the 
diagram may be also calculated and thus the accuracy of the approximation 
method may be judged. 


= 851.5 kg s/i 


m 


(8) 
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Thus is obtained for full weight cancellation (X = l) from the 
formulas of my prior report (ref. 5) the following 


fg = 0.1635 m'l 
f* = 0.1380 m I 
Pg = 5960 kg 


(9) 


The approximation calculation was conducted from f = 0 to f = 120 mm 
in 6 steps of = 20 ram according to the equation 


= -0.08515 + ^0.00725 + v?_i - 0.00020(Pi_i + Fi) 


( 10 ) 


with which numerical table 1 was obtained. 

Since, at f = 120 mm, we have already approached toward the maximum 
force P, the steps following are reduced. From f = 120 mm on, calcu- 
lation is then made with Af = 4 mm according to the equation 


= -0. 01703 + yjo. 00029 + yf_i - 0.ooooJ4-(Pj^_]_ + F^) (11) 

The citation of the individual values is superfluous, since the calcula- 
tion is made basically in the same manner as carried out in the numerical 
table 1. 


For f = 160 mm is obtained v = 0.755 m/s. The step following 
leads to a negative value below the root sign of equation (ll), so that 
we should calculate with the sectional formula (3) with 


A = 882 mkg 
Jn-1 = 853.5 nikg 

^n-1 = 577^^ kg 
Fn.i = 5130 kg ^ 


Then we obtain 


Af. 


n 


^ 882 - 853.5 

5^52 


= 0.0052 m 


(13) 


and, consequently, fg = O.1652 m. Since, according to the discussed 
method, the maximum value of the force is given as Pg = 6012 kg (fig. 3)> 
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the deviation of this from the true value according to equation (9) is 
+52 kg = +0.9 percent, whereas the error in the determination of the 
maximum spring stroke is +1.7 = +1«0 percent. In the example 

treated, the accuracy of the approximation method is therefore excellent, 
especially when considering that for its completion altogether only 
17 steps were calculated for which, inclusive all side -calculations, not 
even one-half working day was required. 

The approximation calculation with quadratic damping may he made 
still faster since in that case the damping curve shows an almost linear 
course. Beginning with the same initial danping as with linear damping, 
calculation is made from f = 0 to f = l60 mm in 8 steps with 

Kp = 202.7 kg s^/m^ according to the formula 


Vi 


0.98033 ^vf.i - 0.0002o(Pi_i + Fi) 


ilh) 


Because the next step to f = I80 mm would lead to a negative value 
below the root sign, from f = 160 mm on we use, with 


A = 882 mkg 

Jn-1 = 781.5 mkg I 

Pn_i = 55*^6 kg 

Fn-1 = 5130 kg ^ 

the end-section formula (3) from which is obtained 


( 1 ^) 


Afn ~ 0.0188 m ( 16) 

Thus we have fg = O.1788 m and Pg = 5676 kg. In exact calculation one 
obtains from the equation (ref. 'J , eq. (l4)) 

tg - 6.7914 - 3 .i 46 o = 0 ( 17) 

by trial and error, i|fg = 5.2629J therefore, 

f = 5.2629 -% = 0.1780 m (18) 

With this we obtain 

Pg = 491 + 29000 X 0.1780 = 5653 kg 


(19) 
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The approximation calculation for quadratic damping with an error of 
+0.8 mm = +0.45 percent in the maximum spring stroke and an error of 
+23 kg = +0.41 percent in the maximum force is, therefore, even more 
exact than for linear damping. 


As a third example for the approximation calculation, a pneumatic 
spring strut with linear damping shall be investigated. Said calculation 
will employ the following basic values 


M = 100 kg s^/m 
Fq = 500 kg 
k = 1000 kg s/m ^ 
Vq = 3.8 m/s 
A = 722 mkg 


(20) 


The equation for the compressive line reads (with 


f in mm) 


F = Fq( 1 - f/200) 


3 


( 21 ) 


Again 6 steps of Af = 20 mm are calculated from f = 0 to f = 120 mm 
according to the equation 


v^ = -0.1 + ^0.01 + vf.i - 0.00020 (Pi_l + Fi] (22) 

and from here on, the succeeding steps with = k mm are calculated 
according to the equation 


Vi = -0.02 + 


y 0.0004 + - o.oooo4(Pi_i + fJ ( 23) 


At fn-1 = 172 mm the residual calculation now sets in which leads to 

Afn 722 - 719.7 — ^ 0.00035 mm s: 0.4 mm (24) 

1(6709 + 6442) 

Thus, we have f_ = 172.4 mm and P_ = 67IO kg. Of the total work 

O 0 

A = 722 mkg, the airstrut absorbs 271 mkg = 37*5 percent while 
451 mkg = 62.5 percent dissipated by the oil damping. This distribution 
markedly differs from that of figure 3; with a straight-line spring 
characteristic the spring absorbs here 53 • 1 percent of the total energy 
and the linear damping 46.9 percent. 
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In elaboration of these investigations, the spring compression 
time tn is also to be determined, which elapses until attainment of 
the maximum loading deflection. 

Generally the following expression applies 


V 


dt 


or dt 


V 


(25) 


from which is had by integration 

fi 



0 


( 26 ) 


Since this expression becomes indetermined for v = 0, that is, for 
^i = fn> its value is found by step-by-step calculation to fn-1 and 
exclusive application of a residual formula. 


To set up said residual formula, the law for the velocity course in 
the last increment must be known or must be assiomed. Here is assumed 
that for linear, as well as quadratic damping, between fn-l f^ 

the velocity decreases per figure 5j this assumption will be more or less 
found to correspond to actuality. Then we have 


1 

V 



(27) 


and, with it 


^i = “tn-l + 



^n-1 + 


2 Af^ 

Vn-1 


2\/Afn(fn - fj) 
vn-1 


For i*! = fn is thus obtained 


( 28 ) 




= tn_i + 


2 Af n 
Vn-1 


(29) 


The step-by-step calculation of the oscillation time, for the example 
per figure 3, with linear damping, yields tjj_]^ = 0 . 06 k 6 s, so that with 

Afn = 0.0052 m and — i— = 1.3238 s/m, we obtain t^ = 0.0781f s. The 


Vn-1 
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oscillation time t^, however, is equal T/4 when T denotes the oscil- 
lation period. This is obtained as 

T = ^ = — 0.3811 s (30) 

whereas 4tj^ = O.3136 s. Thus, the value is by 0.06T5 s = 17.7 percent 
below the true value. 

For quadratic damping an error is found of a similar magnitude. 
According to Klotter (ref. 1, p. 95) the oscillation period at quadratic 
damping is nearly equal to that with complete absence of damping 


T « ^ = 0.3685 s 

03 


(31) 


In the step-by-step calculation of the time intervals 
tn-i = 0.0561 s so that, with Af^ = O.OI88 m and 


, we arrive at 

= 0.6982 

Vn-1 


s/m. 


we obtain t^j = O.O823 s. Thus we have here also 4 tn = O.3292 s which 
is, by 10.7 percent of the correct value, too^ small. 


Although the time determination was made with the greatest care, 
the results are unusable. The time may be calculated with sufficient 
accuracy then only when the velocity as a function of the spring stroke 
is given as a mathematical expression. Arbitrary assiunptions will also 
lead to great errors even when the assxiraptions depart very slightly from 
reality. 


In general, the time is of little interest in drop testing of air- 
plane shock struts so that the failure of calculations may not be decried. 
This question was aired here for the sole purpose to show that further 
activities in this field may be pronounced in advance as futile. 

The above treated examples point, with closer scrutiny, to a second 
method of approximative calculation of drop tests on which shall be 
reported below. Namely, when plotting the dimensionless form v/vq 

as a function of f/fg we obtain (see fig. 6) for the linesir spring 

chart and linear damping curve according to figure 3» in solid lines, and 
for the nonlinear case the broken line curve, shown in figure 4. Both 
cTirves run so close together that an attempt was made to replace them 
with a single curve and to give the latter general validity. The approxi- 
mation curve of the course of the dot-and-dash line follows the law as 
shown below. 
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V = VQsJl - f/fg (32) 

At first we shall investigate under what condition the above law of 
equation (32) is truly fulfilled. The attempt, by series expansion of 
the In and arc tan expression of the exact solution, to interpret 
equation (32) for a linear spring characteristic and linear damping 
Csee ref. 5, p. 130, eq. (l2)) had completely failed. (Translator's 
remeurk: Due to a missing comma, the German sentence is difficult of 

interpretation. We have done our best in translating the involved 
sentence, but recommend caution in accepting it.) When considering only 
the first member in both of the series, one already arrives at an equa- 
tion of the 3rd degree between v and f , with which nothing can be 
undertaken. Irrespective of same, the series development for the arc tan 
expressions is unreliable fro m the st art since, per exanple for a = 0 

ipo v/ 1 - - 

and -9 = 0.25, the value =3*87 is essentially greater than 1. 

^cpo - O’ 


A second attempt for solving equation (32) starts out with the 
differential equation. (See ref. 5, p. 130, eq. (3).) 


+ 2 - 3=0 

dt 9 


(33) 


When now making the assumption which, at first appears quite arbitrary, 
that 


2\J2 ^0 


(34) 


that is, calculating with a direct proportionality between d and q> 
and assuming that the value d = 0.25 is attained at the mean 


ratio 


then this assxamption, although distasteful to the 


mathematical expert, should nevertheless offer no great practical errors 
for medium values of d: 


Equation (33) written with (34) gives 


^ ^ ^ -2- = 0 
dt 9 y/2 9o 


when inserting into this equation 


(35) 


9 = 9q\/i - tAg 


(36) 
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and its differential quotient 


dcp ^0 


cpo" 


\/l - ilf/ig 


2<Ptc 


we thus obtain 


. + t - a 52^ , 0 

2i|rg \/2\<P0/ 


For cp = 0 ve have ^Ir = ^Ir_j therefore, 

& 

- ^2_ + * - a = 0 

2+g 


or 




When inserting equation (39) into ( 38 ) we obtain 


Furthermore, since according to equation (39) 


we obtain finally 


or 




-1 + i/tg + iW%f = 0 


fp = 9o\/i - t/ig 


(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


This is the assumed velocity co\arse per (32) or ( 36 ). With this we again 
return to equation (32). When assuming generally such velocity course 
with the compression of the spring strut, a second method is obtained for 
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the approximative calculation of the force-deflection diagrams. For a 
linear spring characteristic and linesir damping with equation ( 32 ) we 
have 


P = Fq + cf + kvo\/l - f/f| 


g 


(*^5) 


By setting the differential quotient to zero, one obtains 


^ = 0 = c - 
df 2f 


g \/l - f*/fg 


(»^6) 


or 


= fg 1 


{ kvoV 

When inserting this expression into equation (^5), we obtain 


k^VQ^ 

O 


or, in the form written per my former report (ref. 5) 

2 2 
■a 9o 

S = CT + i|;g + 
g ® ^l^g 


(47) 


(48) 


(49) 


The maximum value of the spring deflection fg is determined from the 
energy balance. First we obtain with x = f/fg 




kvdf = kvQfg 
0 0 




X dx = j kvQfg 


(50) 


This implies that the mean velocity v„ between f = 


0 and f = f _ is 

O 


equal vq. With 


Jd = I ^^ 0^8 


( 51 ) 
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is, furthermore 


Jjr + Jj) - 2 ^0^ 


g 


( 52 ) 


When combining the two last equations we obtain 


+ 2f 


I kvQ - (Mg - L - Fq) 


■g -g 

and finally 

I kvQ - (Mg - L - Fq) 


■ (^)' = 


( 53 ) 


2 /I ^"^0 - (Mg - L - Fq)'^ 


"g= ■ 




) 


(51^) 


Written in a different form, equation (5^) reads 


ik = 
^g 


- (I ^'Po - «') + \j^0^ + 


( 55 ) 


When calcixlatlng the approximate values S„ by equation (i<-9) and 

. Cj 

tg by equation (55); and comparing them with the exact values (ref. 5, 
p. I5O; eqs. (15a) and (20)), it is found that with freedom from damping 
(■Q =0) both values show perfect agreement. For •8 = 0.2 and 0.4 fig- 
ures 7 and 8 apply. It may be seen that for a = 0 and a = -1, that 
is, for 1 < A + a % 2, corresponding to the actual conditions with 
CT ss 0.5 and k = 1, the errors of Sg and i|fg remained below 1 percent 
and in the spring stroke determination for a = 0 and 8 = 0.2 only 
rose to 1.6 percent. For a = 1, the errors increase with decreasing q)Q 
which regions, however, are practically meaningless since, in the first 
place the spring struts are always preloaded, and. secondly, cpQ hardly 
ever is smaller than 5 (ref. 5, p. I3I, left colimm). 

Thus, calculations may be made with the second approximation method 
at least as long as more exhaustive test resiilts from experts will place 
the once-for-all necessary preliminary calcvilations of oil-air spring 
struts on fully reliable foimdations. 
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In the next partial report the here obtained results shall be used 
in order to gain a general persi>ective on the maximum forces, maximum 
spring deflections and spring capacities occurring in the landing shocks 
of oil-air spring struts. 


Translated by John Perl 
Lockheed Aircraft Corp. 
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TABLE I 

EXAMPLE OF AN APPROXIMATION CALCULATION 



F 

kg 

V 

m/s 

kiv 

kg 

P 

kg 

J 

mkg 

0 

491 

4.2 

5,576 

4,067 

0 

20 

1,071 

3.992 

5,399 

4,470 

85.4 

4 o 

1,651 

5.751 

3,194 

4,845 

178.5 

6 o 

2,251 

3.473 

2,957 

5,188 

278.8 

80 

2,811 

5.151 

2,685 

5,494 

585.6 

100 

5,391 

2.771 

2,560 

5,751 

498.1 

120 

3,971 

2.311 

1,968 

5,939 

615.0 
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P 



Figure 1.- Cutout of the force -deflection diagram between fj.j^ and 
f| and end-cutout between f^.^ and f^ (linear damping). 


P 



Figure 2.- Cutout of the force -deflection diagram between and 

f- and end -cutout between fri-i % (quadratic damping). 

1 ^ 
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Exact volues to 
maximum force 



Figure 3.- Example of the step-by-step calculation of the force -distance 
diagrams with linear and quadratic damping. 
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0 20 40 60 80 100 120 140 160 172.4 

f mm 


Figiire 4.- Example for the step-by-step calculation of a force -deflection 
diagram with an undetermined spring chart and velocity -proportional 
damping. Oil -air strut. 
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Figure 5.- Approximation calculation of the last time -period. 



f/fg 


Figure 6.- Course of the spring compression velocity in function with 

the spring stroke. 
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^0 


Figure 8.- Exactness of the approximation formula (equ. 49) for determina- 
tion of Sg. 
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THIRD PARTIAL REPORT - THE LANDING IMPACT 
OF OLEO-PNEUMATIC SHOCK ABSORBERS* 

By K. Schlaefke 


Siinnnary: The stress-stroke diagram of an oleo -pneumatic leg under 

the first impact at landing can he obtained by superposing onto the dia- 
gram of the undamped compressed-air leg, the damping diagram according 
to the second approximate method described in the preceding partial 
report (ref. l) . The maximum stroke, maximum load, and efficiency to be 
expected for various damping, lift coefficients, and heights of drop are 
computed and represented diagrammatically for practical use. Lastly, 
it is indicated how the newly introduced damping factor for oleo- 
pneiunatic legs can be determined by experiment. 

The present report may be joined to the second partial report 
(ref. l) without further introduction. First to be considered is the 
compressed-air strut without damping, a problem which, in view of the 
generally known principles of thermodynamics, presents no difficulties. 

The natural frequency co has proved itself useful as characteristic 
value of linear spring systems because it enables the investigations to be 
represented in generally applicable dimensionless form. But, since <u 
is not a constant for nonlinear systems, a different characteristic 
value must be looked for. As such, the theoretical maximum spring 
travel or stroke fQ appears to be the most suitable for oleo-pneumatic 
struts. It is the theoretical stroke required to raise the force F from 
the initial tension Fq to F = w. On actually constructed struts fQ 
ranges from 230 to 600 millimeters if one assumes, besides the limits 
indicated by Michael (ref. 2), ® 0.75 for the maximum travel fg. 

The compression curve follows the law 

F(fQ - f)^ = FQfQ^ (1) 


*"Zur Kenntnis der Kraftwegdiagramme von Flugzeugfederbeinen, 

3. Teilbericht: Der Landestoss von Slluftfederbeinen." Technische 

Berichte, Bd. 11, Heft. 5, May 15, 19*^^^, pp. 137-1^1. 
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with 





( 2 ) 


equation (l) becomes 


Sp = a(l - ilO“^ = 


(3) 


The energy absorbed by the compressed-air spring is 



Spdilf = a 


(1 - . 1 

k - 1 


o<l>2 




the rise of the compression curve 


— = ak(l - iir)“(l^+l) 
d>lf 



(5) 


The variation of the polytropic curve (fig. 1) should not be much 
different from that of the adiabatic curve because the landing impact 
is generally too fast for any appreciable heat flow. In the following, 
k = 1.3 is assumed, which incidentally is the same figure used by Inner 
(ref. 3) and Hadekel (ref. 4) on the basis of experiments. 

With k = 1.3> table 1 and figure 2 are obtained for the functions 
Ol, ^ 2 > equations (3) to (5) with respect to It should 

be noted that the expansion of, for example, in series 


= (1 - t)**^*^ = 1 + 1.3i + l.i*-95^2 + i.6445ilr3 + 


( 6 ) 
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is useless because this series converges very slowly for practical 
values of t . For \[f = 0.5 the error of the series value when breaking 
off after the term of third degree is already almost 10 percent, while 
for =0.8 the value still amounts to less than half of the correct 
value . 

The compression fj. under static load is obtained from the equation 


Mg(fo - = Fofo^*^ (7) 

at 

fj./fo = 'kr = 1 - (8) 

The relationship between and a is shown in figure 3* 

According to Hadekel (ref. U) fg/fQ mostly equal to 0.75 and 
fr/fg = "to 0.667; hence ^Ifj, = 0.375 to 0.5 and a = O.Uo to 0.5^ 

according to equation (8) or figure 3* In the following a = 0.5 is 
used. This is in full agreement with Michael (ref. 2) who puts 
Spg = 2.2 to 3*5^ while table 1 gives Spg = 3 *03 for a = 0.5 and 

fg/fo = 0.75. 

The energy balance for the landing impact of the undamped compressed- 
air leg (6 = O) reads 



F df = Mg(H + (1 - k)fg) 


(9) 


or, in dimensionless form 



Sp dt = ^ + (1 - X)tg 


From the last equation follows 



(l - X.)\|fg 


( 10 ) 


8=0 


a 


0.3 


( 11 ) 
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Equation (ll) can be solved only for X = 1 with respect to ^ g. 
Therefore, it is better to simply consider ''l^g as independent and 
n/fo as dependent variable (table 1, columns 5 and 6). For round 
values of H/fo the corresponding values of 'Ifg are obtained by 

interpolation. Incidentally, it should be noted that the maximum 
travel for the undamped spring impact (X = O) and H = 0 is 
\(fg = 0.68066, a phenomenon already explained elsewhere (ref. 5)* 

The efficiency t] for the pure compressed-air leg is 

(1 - - 1 _ (1 - Og) ■ (i - 

' " 0.3.,(l - ■ 

by which equation column 7 of table 1 was computed. 

This concludes the study of the undamped compressed-air leg. 


TABLE 1.- FUNCTIONS FOR COMPUTING THE COMPRESSION CURVE 


1 

2 

3 

4 

5 

6 

7 


$1 

0 0 

$0 

(H/fo) 

(H/fo) 






0 

II 

0 

II 

cO 

6 = 0 , X = 1 

n 

0 

' 1 

0 

1.3 

0 

0 

1 

0.1 

1 . 1466 , 

0.1070 

1 . 6564 

- 0.0465 

0.0535 

0.933 

0.2 

1.3365 

0.2306 

2.1719 

- 0.0847 

0.1153 

0.863 

0.3 

1.5899 

0.3764 

2.9527 

-0.1118 

0 .1882 

0.789 

0.4 

1.9427 

0.5520 

4.2091 

- 0 . 124 o 

0.2760 

0.711 

0,5 

2.4623 

0.7704 

6.4020 

- 0 . 11^48 

0.3852 

0.626 

0.55 

2.8237 

0.9024 

8.1575 

-0.0988 

0.4512 

0.582 

0.6 

3.2910 

1.0546 

10.696 

-0.0727 

0.5273 

0.534 

0.65 

3. 9149 ' 

1.2340 

i 4 . 54 i 

-0.0330 

0.6170 

0.485 

0.7 

4.7835 

1.4500 

20.729 

0.0250 

0.7250 

0.433 

0.725 

5.3563 

1.5766 

25.321 

0.0633 

0.7883 

0 . 4 o 6 

0.75 

6.0629 

1.7190 

31.528 

0.1095 

0.8595 

0.378 

0.775 

6.9528 

1.8814 

40.173 

0.1657 

0.9407 

0.349 

0.8 

8.1033 

2.0690 

52.671 

0.2345 

1.0345 

0.320 

0.825 

9.6394 

2.28961 

71.607 

0.3198 

1.1448 

0.288 

0.85 

11.7783 

2.5556 

102.08 

0.4278 

1.2778 

0.255 

0.875 

14.9285 

2.8670 

155.26 

0.5585 

1.4335 

0.219 

0.9 

19.9526 

3.3176 

259.39 

0.7588 

1.6588 

0.185 

0.925 

29.0015 

3.9170 

502.70 

1.0335 

1.9585 

0.146 

0.95 

49.1291 

4.8550 

1277.4 

1.4775 

2.4275 

0.104 

0.975 

120.970 

6.7476 

6290.4 

2.3988 

3.3738 

0.057 

1 

00 

1 °° 

00 

00 j 

00 

0 
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The force-deflection diagram of the oleo shock ahsorher is 
obtained by superposing the damping diagram on the static diagreun 
of the pneumatic shock absorber. 

The damping force D proportional to the velocity follows in 
good approximation the law 


D = Dq\/i - f/fg 


(13) 


as proved in the previous partial report (ref. 1, equation 32 , and 
fig. 6). 

With a view to employing the characteristic value fQ for charac- 
terizing the damping also, is defined as the damping force that 

occurs in a landing impact from a drop of height fQ 


DfO = k sjigfo 

(14) 

Furthermore, since 


Dq = kV^ 

(15) 

the elimination of k from equations (l4) and (15) 
equation (13) leaves 

together with 

D = DfQ^/H/fQ \/l - f/fg 

( 16 ) 


or, with 



(17) 


in dimensionless form 


g 


( 18 ) 
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According to the second partial report (ref. 1, equation (50)) 


Ad = 



Sp d* = I 5^tg 


(19) 


with which the energy balance that serves for defining becomes 




Q- - * 6 ) 


-0.3 


- 1 2 


0.3 


+ I = ^ + (1 - \)ig (20) 


For the numerical calculation, equation (20) is best solved with 
respect to H/f^ as unknown. By equation (ll) 


H 


3 


+ 



( 21 ) 


Tlie root disappears for the values of 6 and ’I'j, given in 
table 2. 


TABLE 2 

EVALUATION OF EQUATION (21) 


6 

Vg 

H/fo 

Figure 4 

1.5 

0.5805 

0.0842 

Point A 

3 

0.3456 

0.1194 

Point B 


The curves have horizontal tangents in the two points A and B 
(fig. 4). But it is to be assumed that here the approximate calculation 
does not agree with reality; at any rate, there is no physical explana- 
tion for this variation of the curve. However, since such small values 
of H/fo are practically of no significance, the bare hint at the limits 
of the approximate calculation should be sufficient. 
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Thus, equation (21) is represented in figure 4 with the reservation 
that the anticipated maxinnam compression can be read off in terms of 
6, X, and H/f(). 

In figure the maximum compression for free-drop and weight- 
balanced impact are compared; then Q( ig) = =l)/^g(X=0)‘ Figure 5 

similarly serves for defining the compression with weight -balanced impact 
from the experimental values of the free-drop impact, as explained in 
reference 5 for linear spring characteristic curves. 

Before proceeding to the calculation of the maximum force, the 
limits of H/fQ and 6 up to which the calculation must be carried 

out so as to include the partial range, are determined. 

It has been established that Tq ranges between 230 and 600 milli- 
meters . In the drop test the maximum height of drop H is about 1000 
millimeters, equivalent to a maximum rate of drop Vq = 4.43 meters per 
second. Thus H/Tq is at the most equal to 1.67 for large aircraft and 
rises to a little more than 4 for small airplanes. Consequently, it is 
sufficient to extend the calculation to H/fo = 4, as already done in 

the evaluation of equation (21) in figure 4. 

To define the average damping factor 6 to be expected, recourse 
is had to the damping magnitude for linear spring characteristic curves, 
which as a rule lie at d = 0 . 2 ^ (ref. 6). 

Equating the two initial dampings, the equation 

2^<P0 = = 2df/fg^ = 6yfe/fQ (22) 


with n = 2.2 to 3.5 (ref. 2, p. 395) and fg/fo = 0.75 gives 

0.5/474 o.5\/? 

5 = to = 1.2 to 1.5 

Vo?75 


(23) 


Accordingly, the calculation was made with S = 0, & = 1.5, and 

6 = 3^ which surely includes the practical flight range. 

The maximum force S_ is obtained by differentiating the equation 

o 

S = Sp + Sp = a(l - t)’ + 6/H/fo/l - ^/i'g (24) 
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with respect to t and putting the differential quotient equal to zero. 
Hence 


(1 - 


This equation must be solved by trial. It is best to assume 


as independent and 6 


t/’ 


S 


as dependent variable, and ultimately define the 
g 


applicable value of the correct value 6 by interpolation 


For H/fo = 1 . 5 ^ X = 1 , a = 0.5, and i|fg = 0 . 6 t 6 , a table such as 
table 3 is obtained. Interpolation for 6 = 1.5 results in 


^*/^g = 0.9666 
= 0.6536 ‘ 
Sg = 2.319^ 


(26) 


TABLE 3 

EXAMPLE FOR SOLVING EQUATION ( 25 ) 


II 

0.96 

0.97 

0.98 

5 = 

1.594 

1.444 

1.234 


Figure 6 represents the drop diagrams for several values of H/fo 
for 6 = 1.5 and complete lift relaxation, that is X = 1 . The point 
defined by equation ( 26 ) is indicated by the letter A. It is seen that 
the maximum value formation starts expressly between H/fo =1 and 
n/fo = 1 . 5 , but that this maximum exceeds the initial force only when 
n/fQ still has risen a little above 1 . 5 * 

This is clearer yet in figure 7 . The parabolas starting in the 
point with the coordinates H/fo = 0 , Sg = O.5 represent the initial 
forces. For S = I.5 and X = 1 the steeply rising branch is valid 



NACA TM 1373 


• III-9 


from H/fo = 1.53 on. For such marked damping as 6 = 3 the maximum 
force for the weight-balanced impact in the entire range in question is 
equal to the initial force, while for the undamped impact from H/fo ^ 3 
on, it follows the steep curve. 

In figure 8 the maximum force with weight-balanced and free-drop 
impact is compared. In contrast to the conditions for straight spring 
characteristic curves, the quotients Q(Sg^ = Sg( with about 

0.4 are now only half as great as for the linear curve (ref. 5, fig. 7). 
This phenomenon is comprehensible without further explanation by a glance 
at figure 6. 

The efficiency q follows from the equation 


the interpretation of which is given in figures 9 and 10. The damped 
diagrams in the practical range of damping are much fuller than the 
undamped ones. Thus, for H/fQ = 2 and 6 = 1.5 the ratio 
Q(q) = 0) “ explained by figure 6: for 

H/fQ = 1.5 the efficiency with 0.95 is nearly equal to unity, whereas 
the diagram is almost a rectangle. At increasing values of i^g, q 


decreases rapidly; but to the free-drop impact there always corresponds 
a greater spring travel or stroke tg than to weight -balanced impact 
(figs. 4 and 5). 


The energy absorbed by the strut at the first landing impact is 


A = H/fQ + (1 - X)tg 


(28) 


whence by equation (19) 



^\/H/fotg 

H/fo + (1 - x)ig 


(29) 


As indicated by the representation of equation (29) in figure 11, 
at 6 = 1.5 air and damper oil participate about equally in the energy 
absorption, and which is fairly independent of the height of drop H in 
the practical range. 
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Figure 11 and equation (29) show the way by which the damping 
factor 5 can be determined by a test. 

With full lift relaxation 


s/Wo 



(30) 


Simply measure the maximum compression for a given height of drop 

O 

and calculate Ap by equation (4) or table 1. Then A can be 
computed by equation (28) and with Ap = A - Aj. the damping factor 6 
obtained by equation (30). This calculation is further facilitated by 
figure 12. 

With it, the present report is concluded. Whether I am successful 
in extending the described approximate method to the calculation of oleo- 
pneumatic legs with tires, I don’t know. However, the present three-part 
total report should, I hope, give not only hints to landing -gear diagrams 
for construction and test, but also serve to stimulate further research. 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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Figure 1. 


Figure 


Figure 3.- 


GO 



The compression curves of compressed-air legs. 



i.- Fimctions for calculating the compression curve. 



a 


Initial tension and spring travel under static load. 
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Figures 4 and 5.- Maximum compression of oleo -pneumatic legs at landing 

impact. 



Figure 6.- Examples of force -stroke diagrams for b = 1.5 and = 1 
(complete lift relaxation) at several heights of drop H. 
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H/fo 

Figure 11.- Proportion of oil damping energy to total energy. 



I I I I I I 

0 0.2 0.4 0.6 OjB KO 


Figure 12.- Diagram for determining the damping B from the test values 

H, tg and Aj)/A. 
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